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ABSTRACT 

We have used 605 days of photometric data from the Kepler spacecraft to study KIC 6614501, 
a close binary system with an orbital period of 0.15749747(25) days (3.779939 hours), that 
consists of a low-mass subdwarf B (sdB) star and a white dwarf. As seen in many other 
similar systems, the gravitational field of the white dwarf produces an ellipsoidal deformation 
of the sdB which appears in the light curve as a modulation at two times the orbital frequency. 
The ellipsoidal deformation of the sdB implies that the system has a maximum inclination 
of ^40 degrees, with i ps 20° being the most likely. The orbital radial velocity of the sdB 
star is high enough to produce a Doppler beaming effect with an amplitude of 432 ± 5 ppm, 
clearly visible in the folded light curve. The photometric amplitude that we obtain, K\ — 
85.8 km/s, is ~12 per cent less than the spectroscopic RV amplitude of 97.2 ± 2.0 km/s. 
The discrepancy is due to the photometric contamination from a close object at about 5 arcsec 
North West of KIC 6614501, which is difficult to remove. The atmospheric parameters of 
the sdB star, T cS = 23 700 ± 500 K and log 5 = 5.70 ± 0.10, impl y that it is a rare o bject 
below the Extreme Horizontal Branch (EHB), similar to HD 1881 12 (Heber et al.ll2003l) . The 
comparison with different evolutionary tracks suggests a mass between ~0. 18 and ^0.25 M©, 
too low to sustain core helium burning. If the mass was close to 0.18-0.19 M©, the star could 
be already on the final He-core WD cooling track. A higher mass, up to ^0.25 M , would 
be compatible with a He-core WD progenitor undergoing a cooling phase in a H-shell flash 
loop. A third possibility, with a mass between ^0.32 and ^0.40 M Q , can not be excluded and 
would imply that the sdB is a "normal" (but with an unusually low mass) EHB star burning 
He in its core. In all these different scenarios the system is expected to merge in less than 3.1 
Gyr due to gravitational wave radiation. 
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1 INTRODUCTION 

Hot subdwarf stars (sdBs and sdOs) are found in all galactic 
stellar populations and they are the main source of the UV- 
upturn phenomenon in early-type galaxies (Greggio & Renzinil 
1990; Brown et al, I I 19971) . Subdwarf B stars in particular are post 
red giant branch (RGB) stars with thin (<0.01 Mq) inert hydrogen 
envelopes. Most of them have masses close to 0.5 Mq, with a peak 
near 0.47, and are core helium-burning objects on the extreme hor- 
izontal branch (EHB). If the mass of an sdB star is below ~ 0.3 
Mq, the He-core ignition never took place. These rare stars are the 
progenitors of He-core w hite dwarfs. A recent extensive review on 
sdB/sdO stars is given bv lHeberldioO^) . 

In order to reach the high temperatures and surface gravities 
typical of sdB stars, their progenitors must have lost almost the en- 
tire hydrogen envelope near the tip of the red giant branch (RGB). 
This may happen in a close binary after a common-envelope (CE) 
phase, when the orbital angular momentum is transferred to the 
red giant envelope spinning up the envelope until enough energy 
has been accumulated to eject it. Indeed, about half of sdBs re- 
side in close binary sy stems with white dwarfs (WDs) or low-mas s 
main-sequence stars dMaxted et al.ll20oH ; iNapiwotzki et al.ll2004l) . 
The various formation channels of sdB stars have been studied by 
lHanetaljf2002ll2003h . 

Presently, more tha n one hundred close b inaries with an sdB 
comp onent are known dGeier et all l2011al lbl ICopperwheat et all 
and an important fraction must be composed by sdB+WD 
binaries. These systems are potential progenitors of type la super- 
novae (SN), provided that the total mass exceeds the Chandrasekhar 
limit and that the initial separation is small enough to merge in 
a Hubble time due to gravitational wave radiation. The first can- 
didate found to be a type la SN pro genitor is KPD 1930+2752 
dMaxted et afll2000l : iGeier et alj|2007l) , with a total mass of 1.47 
Mq when we assume a canonical value of 0.5 Mq for the sdB com- 
ponent. Another interesting sdB+WD system, KPD 1946+4340, 
was observed by the Kepler spacecraft and studied in great detail 
dBloemen etalfeoi lh . Thanks to its brithness (V=14.28) and high 
inclination (87.1°), KPD 1946+4340 shows very clean primary and 
secondary eclipses and a very accurate light curve modelling was 
possible, including not only ellipsoidal deformation and Doppler 
beaming, but also WD reflection and gravitational lensing. 

The target discussed in this article, KIC 6614501 (alias 
2MASS J19365001+4201436), is another such sdB+WD system. 
It was selected among a list of sdB pulsator candidates inside 
the Kepler field of view, with the main goal of performing de- 
tailed asteroseismic studies o f these stars. During the survey phase 
(Oste nsen et alj|2010bl.l201 ll). it was observed by the Kepler space 
telescope dBorucki et al.ll2010h i n the framework of th e Kepler As- 
teroseismic Science Consortium dGillilandetal.l2010|). The first re- 
sults a fter 1 month of observation are reported by Ipstensen et alj 
d201 ll) : no pulsations were detected and the atmospheric parame- 
ters were determined from spectroscopy at the 4.2 m William Her- 
schel Telescope. A low frequency modulation plus its 1st harmonic 
in the Kepler light curve was interpreted as a binary signature, sug- 
gesting a sdB+WD sy stem because no red e xcess was seen in the 
spectrum of the star. 0ste nsen et al] d201 ll) also noted the pecu- 
liar position of KIC 6614501 in the T c ff/log<7 plane, below the 
EHB, implying a low-mass post-RGB object. A mass of about 0.24 
Mq was suggest ed by c omp aring its position with the evolutionary 
tracks of iDriebe et all dl998l) . implying that KIC 6614501 could 
evolve into a He-core white dwarf. 

This particular configuration may occur when the envelope of 



the red giant is ejected before the core has attained sufficient mass 
to ignite helium. In this case, the RG core may cross the region 
of the EHB stars near its low temperature - low gravity edge, be- 
fore evolving to the WD cooling track. These post-RGB low-mass 
objects are rare and only recently their number has increased signif- 
icantly thanks in pa rticular to the resul ts of the ELM=Extreme Low 
Mass (WD) survey dBrown et al.l20ld) . the Kepler mi ssion, and the 
WAS P (Wide Angle Search for Planets) survey (e.g. iMaxted et al.1 
[201 ]]). The ELM survey has found many new double degener- 
ate systems, increasing t he number of thos e previously discov- 
ered by the SPY survey dKoester et al.ll2009l) . A spectacular case 
is J1741+6526, composed by a ELM WD primary of 0.16 M and 
an unseen companion with a minimum mass of 1.1 Mq: with an 
orbital period of 1.47 hours, it should merge in less than 170 Myr. 
Assu ming an inclinatio n of 60 degrees, the total mass would be 1.7 
M dBrown et alj2012h . 

However, most ELM WDs are concentrated at effective tem- 
peratures below 15 000 K. At higher temperatures, the number of 
known ELM WDs (or ELM WD precursors) is very small, with 
only a few candidates. KIC 6614501 is the hottest candidate in a 
region of the T c g /log g plane almost completely empty and crucial 
to understand the evolution of the stars that have a mass close to the 
limit for He-core burning. The system most similar to KIC 6614501 
is HD 188112, with a ~0.24 Mq sdB primary and a massive 
(>0.73 Mq) unseen companion dHeber et al.l2003l) . If the compan- 
ion was a C/O WD with a mass larger than 0.9 M , HD 188112 
would be the precursor of a subluminous SN la. 

In the next sections we will perform a detailed study of the 
properties of KIC 6614501, based on Kepler data plus spectro- 
scopic follow-up. In the last section of this article, our results are 
placed in the context of the ELM white dwarfs. 



2 KEPLER LIGHT CURVE AND EPHEMERIS 

KIC 6614501 was observed by Kepler in short cadence (SC, 59 s 
sampling) during the following monthly runs ("Q" stands for quar- 
ter): Q3.3, Q5.1, Q5.2, Q5.3, Q6.1, Q6.2, Q6.3, Q8.1, Q8.2, Q8.3, 
Q9.1, Q9.2, Q9.3, Q10.1, Q10.2 and Q10.3 for a total duration of 
458.9 days. Moreover it was observed in long cadence (LC, 29.4 
min sampling) during Q3.1, Q3.2 and Q7 (146.4 days in total). 
The data were downloaded from the Kepler Asteroseismic Sci- 
ence Operations Center (KASOC) websiteQ The data files contain 
Barycentric Julian Dates (BJD), raw fluxes with errors, and cor- 
rected fluxes with errors. The corrected fluxes are computed by es- 
timating and subtracting the contamination by close objects, based 
on their spectral energy distribution. However, the colour-based at- 
mospheric parameters (T e g and log g) of the stars in the Kepler 
database are not accurate for stars with peculiar properties such as 
very hot stars. For this reason we decided to use raw fluxes. This 
choice may affect the relative amplitude of the orbital modulation 
and its 1st harmonic: a discussion of this aspect is given in section 
4.1. 

For each run the fluxes were corrected for long-term trends us- 
ing a cubic spline interpolation or a straight line in a few cases. We 
removed outliers (>4 x standard deviation), and the fluxes were 
divided by their mean value in order to normalize all runs to an 
average value of 1. 



1 http://kasoc.phys.au.dk/kasoc/ 
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3 SPECTROSCOPY 

On the nights of 12 and 15 August 2010 we obtained 22 spectra of 
KIC 6614501 with the Kitt Peak National Observatory (KPNO) 4- 
m Mayall telescope and its RC-Spec/F3KB spectrograph. The spec- 
tra cover the H^-FL, region with a dispersion of 0.45 A/pix and a 
resolution of 1.6 A. The exposure time was 600 seconds, yielding 
a median S/N=38 per pixel. 

Similarly, on the night of 27 June 2010, we obtained 10 
spectra at the 2.5-m Isaac Newton Telescope (INT) with the 
IDS235+EEV10 spectrograph, with a resolution of 3.0 A, expo- 
sure time of 1600 seconds, and median S/N=42 when resampled to 
the KPNO 4m dispersion. 

The data were homogeneously reduced and analysed. Stan- 
dard reduction steps within IRAF0 were: bias subtraction, removal 
of pixel-to-pixel sensitivity variations, optimal spectral extraction, 
and wavelength calibration based on arc lamp spectra. The target 
spectra and mid-exposure times were shifted to the barycentric 
frame of the solar system. The spectra were normalised to place 
the continuum at unity by comparing with a model spectrum for a 
star with similar physical parameters as we find for the target (see 
section 3.2). 



Figure 1. Top: the Kepler raw light curve after normalisation of each seg- 
ment (upper panel) and after cubic spline (or straight line) correction (mid- 
dle panel). Bottom: amplitude spectrum at low frequency showing the or- 
bital period and its first harmonic. Only in the two upper panels, in order to 
reduce the size of the figures, we used a light curve with a reduced number 
of points by binning to 29.4 min all the short cadence data (while in the 
lower panel and in all our analysis the data were used with their original 
sampling time). 



A first Fourier analysis up to the Nyquist frequency was per- 
formed using only the short cadence data to check for any possible 
presence of high-frequency oscillations. Excluding some known ar- 
tifacts, only a few peaks were found with an amplitude just beyond 
the threshold which we fixed at 4 times the mean value of the am- 
plitude spectrum (4a). The highest peak is at 579.319 d _1 (670.508 
/iHz), with an amplitude of 28.9 ppm, corresponding to 4.3<r. With 
more data from Kepler it will be possible in the future to confirm 
or not the presence of high-frequency oscillations in KIC 6614501. 

The low-frequency analysis was done using all the 676 744 
data points (SC + LC), keeping their original sampling time. We 
verified that the smearing effect introduced by LC data was smaller 
than the amplitude uncertainty. At low frequency, the Fourier am- 
plitude spectrum shows two peaks that correspond to the orbital pe- 
riod and its first harmonic. The latter is the result of the ellipsoidal 
deformation of the sdB star due to the gravitational field produced 
by the companion. Less obvious is the correct interpretation of the 
orbital modulation, which requires also the phase information (see 
section 4). The light curve and the amplitude spectrum at low fre- 
quency are shown in Fig. 1. 

From fitting the orbital period to the Kepler light curve, we 
obtain the following ephemeris: 

BJDtdb = 2455093.33296(31) + 0.15749747(25) E 

which gives the times when the sdB is closest to the Earth (i.e. 
phase in Fig. 4). The reference epoch corresponds to the first time 
at phase of the Kepler data. 



3.1 Spectroscopic radial velocities 

Radial velocities (RVs) were derived with the FXCOR package in 
IRAF using the IL3, H 7 , and H,, lines. All the RV measurements 
are reported in Table 1. Fits assuming a sinusoidal orbital velocity 
curve confirm the orbital period as 0.157497(5) d, fully consistent 
with the photometric orbital period found from the Kepler data. The 
RV amplitude that we find is 97.2 ± 2.0 km/s. Due to the relatively 
long exposure times of the INT data, the measured velocities are 
smeared-out to result in an RV amplitude that is lower than the fit 
by about 2 km/s. For the KPNO data the smearing effects due to 
the exposure time are negligible. When fitting the data of each ob- 
servatory with the period and phase fixed from the all-data fit, we 
found that the RV amplitude is consistent. The system velocity is 
-9.9(2.2) km/s for the KPNO data, and 0.5(3.2) km/s for the INT 
data. This zero point offset also shows up in the position of the in- 
terstellar Call 3993 A line in the mean spectrum of each telescope. 
Although the measurements of the position of this line are not ac- 
curate enough to calibrate this offset, it is clear that the offset is 
site dependent and does not imply the presence of a third body. Our 
best fit is shown in Fig. 2 and the fit parameters are listed in Table 2. 
The orbital period was fixed to the value determined from Kepler 
photometry. 

After the orbital fits, the spectra were shifted to remove the 
orbital motion, before being co-added to obtain high S/N spectra 
(S/N~150) with minimal orbital line broadening, for both observa- 
tories. These spectra were used to derive the atmospheric parame- 
ters of the star. 



3.2 Improved atmospheric parameters 

The high S/N mean spectra from the INT and KPNO, shown in 
Fig. 3, were fitted to a g rid of LTE spectra computed ad hoc for this 
star dHeber et alj 2000). Since there is no detectable helium lines 

2 Image Reduction and Analysis Facility, written and supported by the 
National Optical Astronomy Observatories (NOAO) in Tucson, Arizona 
I http://iraf.noao.edu/ \. 
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Table 1. RV measurements. 



BJD 


v (km/s) 


error (km/s) 


2455375.481497 


13.5 


25.0 


2455375.500704 


-66.8 


15.9 


2455375.520005 


-94.0 


6.6 


2455375.541801 


-48.6 


11.7 


2455375.561105 


-14.8 


10.4 


2455375.579773 


60.7 


5.3 


2455375.601982 


76.1 


8.1 


2455375.624646 


75.2 


5.7 


2455375.649749 


-33.6 


8.0 


2455375.675079 


-102.8 


8.0 


2455420.659287 


71.5 


11.6 


2455420.671652 


35.3 


11.3 


2455420.678830 


22.7 


10.1 


2455420.685997 


2.8 


6.6 


2455420.693178 


-45.0 


7.5 


2455420.705115 


-70.7 


6.8 


2455420.713657 


-104.5 
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2455420.720826 


-98.5 


9.4 


2455420.759432 


-29.5 


9.2 
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6.6 
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45.6 


5.5 


2455420.788526 


71.8 


11.9 


2455423.635478 


74.3 


7.9 


2455423.642873 


80.7 


10.8 


2455423.650043 


87.8 


9.4 


2455423.657195 


67.2 


5.0 


2455423.664441 


55.8 


8.2 


2455423.716534 


-102.8 


8.7 


2455423.724074 


-99.8 


12.9 


2455423.731241 


-85.5 


30.4 


2455423.738429 


-82.0 


40.5 



4 30.6 420.7 430.8 
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BJD - 2455000 



Figure 2. The radial velocity curve of KIC 6614501 using the orbital period 
as determined from Kepler photometry. See the text for more details. 

in the spectra, we fixed the helium abundance to log (Nhb/Nh) = 
—3.0. The grids were computed at solar metallicity and at metal- 
licities reduced by a factor of 10 or 100 relative to solar composi- 
tion. From an evaluation of the line strength in the models when 
convolved with the instrumental profile, we infer, from the non- 
detection of any metal features in our spectra, that the overall metal- 
licity must be significantly sub-solar, at least log[Z] < —1.5. But 



Table 2. Orbital fit parameters using the orbital period from Kepler pho- 
tometry. 



Offset (km/s) -6.5 ± 1.5 

Amplitude (km/s) 97.2 ± 2.0 

Period (day) 0.15749747 



Phase (BJD-2455000) 375.56366 ± 0.00054 




3800 4000 4200 4400 4600 4800 5000 
Wavelength [A] 
Figure 3. The high S/N mean spectra from the INT and KPNO. 



we know that individual elements may be substantially under- or 
overabundant relative to the solar composition, as is typical for the 
sdB stars due to the competing effects of gravitational settling and 
radiative levitation an d to possible radiation-driven wind (see sec- 
tion 5 of EebeHlooi and references therein). In the end we decided 
to use the grid computed at log[Z] = —2.0 for the final analysis. 
The H/3 line was kept out of the fit, partly due to a problem with the 
background subtraction in the KPNO data, and partly due to con- 
cerns regarding NLTE effects in the core, which are not accounted 
for by our models. For the KPNO data set the H,, line was also kept 
out due to a bad CCD column which translates into a broad artifact 
feature in the mean spectrum after applying the orbit correction. 
Our best atmospheric parameters from INT and KPNO spectra are 
summarised in Table 3. 



Table 3. Best atmospheric parameters from INT and KPNO spectra. 



Data 


T c g 


log(g) 


INT 


24 084[65] 


5.732[11] 


KPNO 


23 332[80] 


5.653[12] 


Adopted 


23 700 ±500 


5.70 ±0.1 
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Figure 4. Phase-folded light curve of KIC 6614501 (data points grouped in 
100 phase intervals) and our best fit, obtained considering two components: 
an orbital modulation that accounts for the beaming effect (black dotted) 
and an ellipsoidal modulation at twice the orbital frequency (green dashed). 
The bottom panel shows the residuals grouped in 100 (blue dots) and 20 (red 
line) phase intervals. Colors are available only in the electronic version. 



4 LIGHT CURVE ANALYSIS 

The light curve of KIC 6614501 can be fit using two components: 
the ellipsoidal deformation of the sdB star and the Doppler beam- 
ing. Our best fit to the data with two sinusoidal components is 
shown in Fig. 4. We obtain ~x^ r!id =\.\l . The amplitudes of the or- 
bital beaming frequency and of the ellipsoidal modulation at two 
times the orbital frequency are 432.4 ± 5.3 and 324.8 ± 5.3 ppm 
respectively. A third component due to the reflection effect by the 
illuminated hemisphere of the companion is too faint to be seen in 
the light curve and can be neglected. This component would have 
the maximum at phase 0.5 and is not seen in the residuals (bottom 
panel of Fig. 4). The lack of an observable reflection effect suggests 
that the companion must be a white dwarf. With such a short orbital 
period, an M-dwarf companion would be easily se en (see e.g. the 
folded light curves of KIC 1868650, KIC 9472174 tOstensen et all 
l2010al lbh. KIC 2991403, KIC 11179 657 iKawaler et al]|20ld ) and 
KIC 7335517 lOstensen et ai]|201lh ). 



4.1 Doppler beaming and photometric RVs 

The asymmetry in the ellipsoid al modulation of KIC 6614501 is 
cause d by Doppler beaming dHills & Paid 1 19741 ; lLoeb & Gaudil 
120031) . 

Doppler be aming, also c alled Doppler boosting, was 
first detected by iMaxted et all j2OO0t> in the sdB+WD binary 
KPD 1930 +2752. Recen t detec tions of Doppler b e aming are re- 
ported by Shp orer et al.l d2010|) and Vennes et al. ( 201 lh . using 
data from the ground, and bv llVIazeh & Faiglerl d2010h , who have 
measured the Doppler beaming caused by a substellar compan- 
ion in CoRoT data. In the context of Kepler, Doppler beaming 



was discussed by lLoeb & Gaudil II20031) and IZucker et al.l J20071) . 
Thanks to Kepler 's sensitivity. Ivan Kerkwiik et alj 120101) have ob- 
tained the first photometric radial velocities (RVs) for the two 
A/B+WD systems KOI 74 and KOI 8 1 . For KOI 74 the photometric 
RVs were confirmed a-posteriori by spectroscopic measurements 
(Bloe men et al.l |2012|) . A good agreement between photometric 
and spectroscopic RVs was found also for th e two sdB+WD bi- 
naries KPD 1946+4 340 iBloemen et aT]|201lh and KIC 11558725 
l lTelting et alj2012t> . 

When a source moves at a radial nonrelativistic speed v r with 
respect to an observer, the observed flux F\ is given by: 



F\ — Fn. 



(1) 



where Foa is the emitted flux and c is the light speed. 
B x = B a +B 3X = 3+(2+dlnF A /dlnA) = 5+dlnF A /dln A is the 
beaming factor jLoeb & Gaudi 2003 1 iBloemen et alj|201 ]]), which 
incorporates two different terms: a geometrical term B g = 3, 
not depending on wavelength, that includes a +1 contribution (en- 
hanced photon arrival rate of an approaching source) and a +2 
contribution (solid angle geometrical aberration). And the Doppler 
shift B s \ that can either increase or decrease the flux of an ap- 
proaching source, depending on its spectral energy distribution 
(SED) and on the instrumental bandpass. In our case, for a hot star, 
the optical flux is on the Rayleigh- Jeans tail of the distribution and 
therefore, when the star is approaching, the blue-shift pushes part 
of the flux outside the visibility range giving a negative contribution 
to the beaming factor. 

When equation (1) is used in broadband photometry, we need 
to compute a weighted mean of the beaming factor that takes into 
account the star's SED and the bandpass, which in turn depends on 
the instrumental response function (ejx), atmospheric transmission 
(e a \) and interstellar medium transmission (e im \): 



(B) = 



J e%\ e a \ eim\ A Fa B\ dX 
J £jA e a \ €i m \ A Fa d\ 



(2) 



In our case ei\ is the response function of the Kepler bandpass, 
e a \ = 1 at any A, and e im \ is the interstellar medium transmission 
corresponding to a reddening E(B — V) = 0.06 ± 0.03 mag (95% 
confidence interval). This value was obtained by collecting multi- 
colour photometry from the Carlsberg Meridian Catalog (SDSS r' 
band, ViZieR catalog 1/304), the Kepler Input Catalog (SDSS g,r,i 
and z bands, ViZieR catal og V/l 33) and the UVEX survey (UVEX 
u, g, r, i and Ha ban ds. | Groot et al1l2009l) . and by following the 
method described in Degro ote et al. ( 201 lb . but applied to a grid of 
sdB spectral energy distributions. 

In order to compute the beaming factor and its uncertainty 
from equation (2), we have used a grid of fully metal lin e-blanketed 
LTE atmosphere models for sdB stars (Heber et al. 200(3) with T e g 
ranging from 21 000 to 25 000 K, log g between 5.1 and 5.9, 
metallicity between solar and 1/100 solar and fixed helium abun- 
dance log(NHo/NH) = —3. The beaming factor is found to be 
(B) — 1.51 ± 0.02 (1.52 without considering reddening). This 
value corresponds to a metallicity of 1/100 solar. The uncertainty 
incorporates any value of the metallicity between 1/100 solar and 
solar (metallicity is the main source of error), and an uncertainty 
of ±500 K and ±0.1 dex on the sdB's effective temperature and 
gravity. 

Once the beaming factor is computed, in order to calculate the 
photometric radial velocity amplitude, equation (1) can be written 
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Figure 5. Left: the field of KIC 6614501 in a real Kepler image with the 
position of the pixel array mask used during the Q7 ran. The pixel size of 
about 3.98 arcsec is not constant over time, differential velocity aberration 
results in the coordinates of the pixel mask changing constantly, typically 
by a few arcsec over each quarter. The three pixels actually used to form the 
optimal aperture, i.e. to collect photons from KIC 6614501, are those with 
a red dot (colors available only in the electronic version). Right: the same 
field in a SuperCosmos Sky Survey 3-color BRI image, which shows the 
target and its close companion at about 5.0 arcsec North West of it. Both 
images are 2' x2'. 



K 



cA B 
(B) 



(3) 



where As is the beaming amplitude. From equation (3) we obtain 
a photometric RV amplitude K\ — 85.8 km/s, 11.7% less than 
the spectroscopic RV amplitude of 97.2 ± 2.0 km/s. The reason 
of this discrepancy is the photometric contamination from a close 
object, at approximately 5.0 arcsec North West of KIC 6614501, 
while the Kepler's pixel size is 3.98 arcsec (see Fig. 5). This flux ex- 
cess dilutes the beaming modulation and decreases its amplitude^. 
The uncertainty on the photometic RV amplitude is dominated by a 
systematic error given by the photometric contamination, which is 
difficult to estimate. Then we have two statistical errors associated 
with the beaming factor (±1.3%) and to the beaming amplitude 
(±1.2%), coming from the fit of the light curve. 



4.2 Amplitude of the ellipsoidal deformation 

From the comparison between photometric and spectroscopic RV 
amplitude, we have seen that only 88.3% of the photons attributed 
to KIC 6614501 actually come from this source. This means that 
not only the beaming amplitude but also the ellipsoidal modulation 
amplitude is reduced by the same amount. Thus we estimate that 
the amplitude of the ellipsoidal deformation, after correction for 
the photometric contamination, is 368.0 ±11.7 ppm. 



5 PROPERTIES OF THE BINARY SYSTEM 

The binary system is described by the mass function, that can be 
expressed as follows: 




M,=0.42 Mo 



: 



. . / P X 1 / 3 {M 1 + M 2 ) 2/3 

sin t = Jfi — — 

\2hG) M 2 



(4) 



M,=0.18 M„ 



20 40 60 80 

inclination [deg] 

Figure 6. Mass of the white dwarf secondary as a function of the system's 
inclination for a wide range of sdB masses, from 0. 18 to 0.42 Mq with 0.06 
Mq steps. The red region on the right is forbidden by the lack of eclipses. 
When we consider also the amplitude of the ellipsoidal deformation of the 
sdB primary, we obtain a family of solutions represented by the blue dots 
and the green or violet 1 a uncertainties (colors available only in the elec- 
tronic version). See text for more details. 



In Fig. 6 the mass of the secondary Q is shown as a function of the 
system's inclination for a wide range of sdB masses. We see that 
the mass of the secondary is always larger than ~0.1 Mq. 

In order to constrain the inclination of the system and the 
mass of the secondary, we can consider also the ellipsoidal defor- 
mation of the primary. The amplitude of the ellipsoidal modula- 
tion, equal to 368.0 ± 11.7 (section 4.1), is roughly proportional 
to rss (Mi/Mi) (Ri/a) 3 sin 2 i, where a is the orbital semi-major 
axis and R \ is the radius of the prim ary. A more detailed treatment 
is given bv lMorris & Naftilanl ( fl993h . 

The radius of the primary can be estimated from 
Ri = (GN'h/g) 1 ' 2 , which gives for example ~0.102 R 
for Mi = 0.19 Mq. Using a gravity darkening coefficient of 
0.479±0.010 (see section 5.1) and adopting 0.27±0.02 for 
the linear limb darkenin g coefficient in the Kepler passband 
dClaret & Bloemenl 1201 fh . we can search for the best values of 
i and M 2 that satisfy equation (3) of iMorris & Naftilanl dl993l) . 
These solutions are shown in Fig. 6 as blue dots. The uncertainties 
in green or violet (68% co nfidence) we re computed by means of 
the simpler equation (6) of Morrisl l 19851 (which, in our case, gives 
results almost identical to IMorris & Naftilanl Il993l always within 
2%), taking into account the observational uncertainties in logg, 
K%, amplitude of the ellipsoidal deformation and orbital period. 
For the mass of the primary we assumed an indetermination of 
±0.02 Mq. Despite the large uncertainties that we see in Fig. 6, 
the best values that we obtain for Ad 2 are in good agreement with 



3 Even using the Kepler pixel data it is not possible to exclude the photons 
coming from the close source in a reliable way because the coordinates of 
the pixel mask change constantly in time, see caption of Fig. 5. 



4 From here on the terms "primary" and "secondary" are used to indicate 
the brightest (sdB) and the faintest (WD) star of the binary, without any 
reference to the initially most massive star. 



© 2012 RAS, MNRAS 000,[TJ|l0] 



Orbital properties of an unusually low-mass sdB star in a close binary with a white dwarf 7 



typical WD masses. More importantly, Fig. 6 suggests a relatively 
low inclination i<40°, with a maximum probability near 20 
degrees. This is true for any primary mass and implies a minimum 
mass of ~0.2 Mq for the secondary. The large uncertainties on 
i and M% in Fig. 6 are roughly proportional to the uncertainty in 
logg, which dominates the error budget. A better determination 
of the surface gravity would be important to further constrain the 
system. 

5.1 Markov chain Monte Carlo simulations 

In order to confirm the analytical results of Fig. 6, and despite 
the level of degeneracy, we made an attempt to constrain the pa- 
rameters of the system by modelling the light curve using the 
LCURVE light curve synthesis code and performing Markov chain 
Mo nte Carlo (MCMC) s imulations ( see ICopperwhea t~et~al] |20 1 d 
and iBloemen et al]|201 ll for details of the code). The gravity and 
limb darkening coefficients for the primary were computed using 
the same sdB atmosphere models as used in the computation of the 
beaming factor (section 4.1), assuming E(B-V)=0.06. The gravity 
darkening coefficient was set to 0.479 ± 0.010, accounting for the 
uncertainty in T c ft and log g, while the limb dar kening coefficient 
was o btained from a 4-parameter law (eq. (5) of IClaret & Bloemea 
l201lh . Then we fixed the RV amplitude of the primary to the spec- 
troscopic value, K\ = 97.2 ± 2 km/s, and we assumed a WD 
radius R2 = 0.015 ± 0.003 R0. Actually this value has a much 
larger uncertainty (and there are no eclipses that constrain the radii) 
but this does not affect the results since the light contribution of 
the white dwarf is only ~0.2%. Three MCMC runs were launched 
with a primary mass of 0.18, 0.24 and 0.36 Mq. Inclination, mass 
ratio and beaming factor were kept as free parameters. The results 
give an average inclination near 20 degrees, in agreement with the 
analytical results. Moreover the formal uncertainties indicate an in- 
clination below 50 degrees at a 95% confidence level for all these 
primary masses. However, these uncertainties are not very reliable 
because of convergence problems that are due to the very wide al- 
lowed range in the inclination and the highly unlinear relation of the 
inclination with other parameters such as the mass ratio. From the 
MCMC simulations we obtained a beaming factor of 1.33 ± 0.03, 
confirming that the photometric contamination must be ~12% in 
order to bring this value to 1.51 ± 0.02 as obtained from synthetic 
spectra. 



5.2 Different scenarios for the sdB evolution 

If the mass of KIC 6614501's primary is below ~0.30 Mq, it must 
be a He-core WD progenitor and we can compare its effective tem- 
perature and surface gravity with He-core WD evolutionary tracks. 
We obtain two possible scenarios: i) for Mi near 0.18-0.19 Mq, 
T e fi and log g are compatible with a He-core white dwarf already 
on the final cooling track (after the possible H-shell flashes). In this 
scenario a higher mass is excluded because the smaller radius due 
to the electron degeneracy would be incompatible with the surface 
gravity and with the amplitude of the ellipsoidal deformation. In- 
deed, looking at lPanei et al. I J2007t) , we see that there is only one 
model that, having already reached the final cooling sequence, can 
fit T e a and log g of KIC 6614501's primary: it is the model with a 
mass of 0.1869 Mq between phase C and D, which corresponds 
to an age of about 180 Myr from the end of the mass transfer. 
The evolution is still relatively fast and the match between theo- 
retical and observed T c n and log g is valid only for about 1 Myr. 



Masses smaller than ~0.18 Mq (e.g. IPanei et al. I d2007l) 's model 
with 0.1604 Mq) are ruled out because the star would never reach 
an effective temperature of 23 700 K. ii) For A/i>0.19 Mq, the 
sdB atmospheric parameters are compatible with a low-mass He- 
core WD progenitor that has not yet reached the final WD cooling 
branch and is evolving along one of the cooling sequences related to 
the H-shell flash episodes. For example T c a and log g are compat- 
ible with the 0.196 M ( fourth cooling branch ) or 0.242 Mq (first 
cooling branch) model of Altha uset al.l l l200ll) . The agreement be- 
tween theoretical and observed T c g and log g is valid for less than 
1 Myr or almost 2 Myr, and the age would be, respectively, 290 or 
only 30 Myr. 

If the mass of KIC 6614501's primary is larger than ~0.32 
Mq, there is another possibility: iii) the sdB is a He-core burning 
object belonging to the extreme horizontal branch with a unusually 
low mass, between ~0.32 and ~0.4 Mq . For such l ow masses the 
theoretical sdB evolutionary tracks dHan et aUl2002l) pass through 
the position of KIC 6614501's primary in the T c ft /log g diagram. To 
explore this possibility we computed a set of zero-age and terminal- 
age EHB (ZAEHB and TAEHB) for sdB masses of 0.33, 0.40 and 
0.47 M , with different envelope masses (0.00001, 0.0001, 0.001 
and 0.005 Mq). These ZAEHB and TAEHB tracks are shown in 
Fig. 7. We see from Fig. 7 that, if the mass of KIC 6614501's pri- 
mary was close to 0.33 Mq, the star would be close to the He-core 
exhaustion while, for a mass near 0.40 Mq, it would be close to 
the ZAEHB and much younger. Such low-mass EHB stars are cre- 
ated from a progenitor with an initial mass between 2 and 2.5 Mq, 
for whi ch the He-core ignition happens under non-degenerate con- 
ditions faan et alJ[2002TlHu et alj2007l : |Prada Moroni & Stranierd 
2009). Usually, sdB stars are considered to be post-He-flash objects 
and thus having a mass near 0.47 Mq. H owever, binary popula- 
tion synthesis models bv lHan et alj d2003l) show that sdB masses 
down to 0.32 Mq can occur in a close binary with a WD compan- 
ion through the so-called second common envelope ejection chan- 
nel. Although these low sdB masses occur less frequently than the 
canonical 0.47 Mq, they should not be neglected. 

With the current data it is difficult to say which of the three 
scenarios proposed, He-core WD, He-core WD progenitor in a H- 
shell flash loop, or low-mass C/O-core WD progenitor, is the cor- 
rect one. We note, however, that the first hypothesis implies a sec- 
ondary mass close to the peak of the WD mass distribution: for 
Mi = 0.19 M Q , we obtain i = 21.2° i 1 !;? and M 2 = 0.57±g:§| 
Mq. An sdB mass greater by a factor of 2 would result in a 
higher luminosity by the same factor. Thus, to know the distance 
of KIC 6614501 would help to discriminate between the different 
possibilities. It is interesting to note that, from one hand, a low- 
mass EHB star would be compatible also with T c n and log g of the 
other two objects close to KIC 6614501 's primary in the T c ff/log g 
diagram (Fig. 7): HD 188112 and SDSS J1625+3632. But on the 
other hand, the known trigonom etric distance of H D 188112 im- 
plies a mass of 0.24±o;J? M faeber et al.ll2003l) . which points 
towards a He-core WD progenitor, although a low-mass EHB of 
~0. 32-0.34 Mq can not be totally excluded. 

In all the three scenarios proposed, KIC 6614501 is expected 
to merge in a time shorter than a Hubble t ime d ue to gravita- 
tional wave radiation. Following PaczynskJ dl967l) , the merging 
time varies between 0.9 (for Mi=0.40 and M 2 =1.01 M ) and 3.1 
Gyr (for Mi =0.18 and M 2 =0.54 Mq). 
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6 KIC 6614501'S PRIMARY IN THE CONTEXT OF 
HE-CORE WHITE DWARFS 

Low-mass white dwarfs have masses below ~0.45 Mq and they 
are generally thought to contain inert helium cores. However, a car- 
bon/o xygen {CIO) core i s not excluded for masses down to ^0.33 
M jPanei et alj l200l IPrada Moroni & Stranierol l2009h . Below 
0.3 Mq, helium burning is ruled out and the white dwarfs must 
have helium cores. These objects are also called extremely low 
mass (ELM) white dwarfs. At least for solar metallicities, ELM 
white dwarfs can form only from binary systems. At higher metal- 
licity it is less clear and the possible detection of single low- 
mass white dwarfs in NGC 6791, one of the oldest (<8 Gyr) and 
most metal-rich ([Fe/H] +0.4) open clusters of our Galaxy, ha s 
prompted a heated debate dOriglia et ai1l2006t iKalirai et alj|2007b . 
The evolution of He-core white dwarf mod els with low- or high- 
meta llicity progenitors h as been studied by ISerenelli et al.l ( 2002) 
an d lAlthaus et ail {2009) respectively. 

With a mass between ~0. 18 and ~0.40 Mq, the primary star 
of KIC 6614501 is a precursor of a low-mass white dwarf. If the 
mass is larger than about 0.32 Mq, it may evolve into a low-mass 
C/O-core white dwarf. For lower masses it will become a He-core 
WD. If the mass is near 0.18-0.19 Mq, it may have already reached 
the final ELM WD cooling sequence. 

The ELM white dwarfs (and their precursors) are a rare class 
of stars. A list of these post-RGB low-mass objects from the liter- 
ature is given in Table 4. Actually the list we found was longer but 
we decided to show only those with a high reliability, for which the 
orbital period was measured through radial velocities and/ or high 
precision photometry from space (Kepler). Moreover, we selected 
only those with logi?<7.0, which co rresponds to exclu de those 
with a mass higher than 0.3 M (see lKilic et ai]|201ld. Fig. 8). 
Some new objects recently found by the ELM survey ferown et all 
|2012|) . which do not meet our criteria, are not in Table 4. Simi- 
larly, we have not included the primary star of SDSS J1257+5428 
(Mars h etalj|201lh because of the large uncertainties, in particu- 
lar in log g. In most cases, the masses reported in Table 4 are ob- 
tained from evolutionary tracks and hen ce are mod e l depe ndent. 
For example, the evolutionary tra cks of IPanei et al.1 ((2007) differ 
from those of Driebe et al. ( 1998), in particular for what concerns 
very low masses. IPanei et aLU2007l) have shown that, when element 
diffusion is included, a dichotomy appears for stellar masses near 
0. 17 Mq : below this limit, the envelope is thick enough for residual 
nuclear burning that slows down the evolution. For most systems in 
Table 4 that have a WD companion, only the minimum mass of the 
secondary is known and thus a more massive neutron star or black 
hole companion can not be excluded. The position of the objects 
of Table 4 in the T e ff/log g plane is shown in Fig. 7. In this fig- 
ure we note that KIC 6614501 is the hottest object, in a region of 
the T e ff /log g plane that is still quite empty. Only other two stars 
reside in the same region: HD 188112, at slightly lower T e fr, and 
SDSS J1625+3632 at almost identical T c g but higher gravity. These 
objects are crucial to study the limit between He-core burning EHB 
stars, that will become low-mass C/O-core white dwarfs, and He- 
core WD progenitors. They can be important also to study the com- 
plex evolutionary path that, through the H-shell flashes, leads to the 
formation of He-core (ELM) white dwarfs. 



5 The difference between ELM WDs and their precursors is subtle. In the 
next lines we will simply use the term ELM WDs for all of them, according 
with recent articles. 




25 20 15 10 



T e „/1000K 

Figure 7. ELM white dwarfs in the T c ff /log g plane. These objects are in 
binary systems with a main sequence star, a white dwarf, or a neutron star 
(black, red, or blue symbols respectively, colors available only in the elec- 
tronic version). The two objects in yellow are in globular clusters. The green 
symbol corresponds to SDSS J1448+1342, that could be either the compo- 
nent of a pole-on binary or a single star. More details on the single objects 
are given in Table 4. The evolutionary tracks are from A lthaus et allfeOOlh . 
with a mass of 0.169, 0.196, 0.242 and 0.292 M . For the models with 
a mass greater than 0.17 Mq, which experience at least one thermonu- 
clear H-shell flash, we show only the first and the last cooling branch (those 
with a longer duration and a higher probability to be populated) for clarity. 
Moreover w e have included als o the final cooling branch of the 0.1868 Mq 
model from Panei et al. (2007), which matches very well the atmospheric 
parameters of KIC 6614501's primary. The numbers close to the triangles 
are the ages in Myrs from the end of the mass transfer. The ZAEHB and 
TAEHB tracks for sdB masses of 0.33, 0.40 and 0.47 Mq are reported and 
correspond to a range of envelope masses between 0.00001 and 0.005 Mq. 
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Table 4. ELM white dwarfs. 
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1V1 \ 1 IV J. (V) 


P >Vrl 
r orb lu 


Eel. 6 


Comp. 


ivi|2' -"-'-(0 


Ref 


V209 ui Cen 


V=16.58 


10866 ± 323 


4.34 ± 0.02 


0.14 


0.834 


2 


reborn WD 


0.94 


1 


WASP J0247-25 


V=12.44 


13400 ±1200 


4.75 ± 0.05 


0.23 


0.668 


2 


A-type 


1.48 


2 


SDSS J1233+1602 


19.83 


10920 ± 160 


5.12 ±0.07 


0.17 


0.151 




WD 


^0.86 


3 


SDSS J1741+6526 


18.27 


9790 ± 240 


5.19 ±0.06 


0.16 


0.061 




WD 


>1.10 


4 


SDSS J21 19-0018 


20.00 


10360 ± 230 


5.36 ± 0.07 


0.17 


0.087 




WD 


^0.75 


3 


SDSS J0917+4638 


18.70 


11850 ± 170 


5.55 ± 0.05 


0.17 


0.316 




WD 


^0.27 


3 


SDSS J01 12+1835 


17.11 


9690 ± 150 


5.63 ± 0.06 


0.16 


0.147 




WD 


^0.62 


4 


KIC 10657664 


13.08 c 


14900 ± 300 d 


5.65 ± 0.04 d 


0.37 d 


3.273 


2 


A-type 


2.5 


5 


HD 188112 


V=10.2 


21500 ± 500 


5.66 ± 0.05 


0.24 


0.607 




WD 


>0.73 


6 


GALLEX J1717 


V=13.7 


14900 ± 200 


5.67 ± 0.05 


0.18 


0.246 


1 


WD 


0.9 


7 


SDSS J08 18+3536 


20.48 


10620 ± 380 


5.69 ± 0.07 


0.17 


0.183 




WD 


^0.26 


3 


KTC 06614501 


15.80 c 


2^,700 + 500 


? 70 + 1 


0.24 


0.157 


() 


WD 


^0.5-1.1 




KOI 1224 e 


14.01 c 


14400 ±1100 


5.72 ± 0.05 


0.20 


2.698 


2 


A/F 


1.59 


9 


NGC6121 V46 


V=18.5 


16200 ± 550 


5.75 ±0.11 


0.19 


0.087 





WD 


^0.26 


10 


KOI 8l f 


11.30 c 


17000 ±1300 


5.78 ±0.13 


0.3 


23.878 


2 


B9-type 


2.7 


11 


SDSS J0152+0749 


18.01 


10840 ± 270 


5.80 ± 0.06 


0.17 


0.323 




WD 


>0.57 


4 


SDSS J0755+4906 


20.09 


13160 ±260 


5.84 ± 0.07 


0.17 


0.063 




WD 


^0.81 


3 


SDSS J1422+4352 


19.79 


12690 ± 130 


5.91 ± 0.07 


0.17 


0.379 




WD 


^0.41 


3 


SDSS J1630+2712 


20.04 


11200 ±350 


5.95 ± 0.07 


0.17 


0.276 




WD 


^0.52 


3 


SDSS J0106-1000 


19.8 


16490 ± 460 


6.01 ± 0.04 


0.17 


0.027 





WD 


0.4 


12 


SDSS J1625+3632 


19.36 


23570 ± 440 


6.12 ±0.03 


0.20 


0.232 




WD 


^0.07 


13 


SDSS J1840+6423 


18.76 


9140 ± 170 


6.16 ±0.06 


0.17 


0.191 




WD 


^0.64 


4 


SDSS J1439+1002 


17.81 


14340 ± 240 


6.20 ± 0.07 


0.18 


0.437 




WD 


^0.46 


3 


SDSS J0849+0445 


19.31 


10290 ± 250 


6.23 ± 0.08 


0.17 


0.079 




WD 


^0.64 


14 


SDSS J1443+1509 


18.58 


8810 ±220 


6.32 ± 0.07 


0.17 


0.191 




WD 


^0.83 


4 


PSR J1012+5307 


V=19.6 


8670 ± 300 


6.34 ± 0.20 


0.16 


0.605 




PSR 


1.6 


15,16 


LP 400-22 


V=17.22 


11170± 90 


6.35 ± 0.05 


0.19 


1.010 




WD 


^0.41 


17,18,19 


PSR J1911— 5958 


V=22.1 


10090 ± 150 


6.44 ± 0.20 


0.18 


0.837 




PSR 


1.34 


20 


SDSS J0822+2753 


18.33 


8880 ± 60 


6.44 ±0.11 


0.17 


0.244 




WD 


^0.76 


14 


KOI 749 


10.85 c 


13000 ±1000 


6.51 ±0.10 


0.23 


5.189 


2 


A 1 -type 


2.2 


11,21 


NLTT 11748 


V=16.7 h 


8690 ± 140 


6.54 ± 0.05 


0.18 


0.235 


2 


WD 


0.76 


22,23,24,25,26 


SDSS J1053+5200 


18.93 


15180 ±600 


6.55 ± 0.09 


0.20 


0.043 




WD 


^0.26 


14,3 


SDSS J1512+2615 


19.24 


12130 ±210 


6.62 ± 0.07 


0.20 


0.600 




WD 


^0.28 


3 


SDSS J0923+3028 


15.63 


18350 ± 290 


6.63 ± 0.05 


0.23 


0.045 




WD 


^0.34 


3 


SDSS J1234-0228 


17.86 


18000 ± 170 1 


6.64 ± 0.03' 


0.23 


0.091 




WD 


^0.09 


13 


SDSS J1518+0658 


17.46 


9810 ± 320 


6.66 ± 0.06 


0.20 


0.609 




WD 


^0.58 


4 


SDSS J1436+5010 


18.23 


16550 ± 260 


6.69 ± 0.07 


0.24 


0.046 




WD 


^0.46 


14 


SDSS J065 1+2844 


19.1 


16400 ± 300 


6.79 ± 0.04 


0.25 


0.009 


2 


WD 


0.55 


27 


PSR J02 18+4232 


V=24.2 


8060 ± 150 


6.9 ± 0.7 


0.2 


2.029 




PSR 


1.6 


28 


SDSS J1448+1342 


19.22 


12580 ±230 


6.91 ± 0.07 


0.25 


i 








3 



Notes: a system's magnitude; b eclipse: 0=not detected, l=secondary eel. detected, 2=secondary and primary eel. detected; c from Kepler Input 
Catalogue; d a secondary solution from the same authors gives T cff /K=14600 ± 300, logg=5.50 ± 0.02, M/M Q =0.26; e KIC 6606653; 
f KIC 8823868; " KIC 6889235; h from USNO; 4 T off /K=17470 ± 750, log g=6.38 ± 0.05 from Liebert et al. 2004; 1 can be either a 
pole-on binary or a single ELM WD, see Brown et al. 2010 for more details. 
References: 1 Kaluznv et al. 2007 , 2 Maxtedetal. 2011 . 3 Brown et al. 2010 . 4 Brown et al. 2012. 5 Carteret al. 2011. 6 Heber et al. 2003 , 

7 Vennesetal. 2011. 8 this paper, 9 Breton etal. 2012. 10 O'Toole et al. 2006. 11 van Kerkwiik et al. 2010. 12 Kilic et al. 2011b. 
13 Kilic et al. 2011a. 14 Kilic et al. 2010a . 15 van Kerkwiik et al. 1996. 16 Callanan et al. 1998. 17 Vennes et al. 2009. 18 Kilic etal. 2009 . 
19 Kawka et al. 2006 , 20 Bassa et al. 2006, 21 Bloemen et al. 2012. 22 Kawka & Vermes 2009. 23 Kawka et al. 2010 . 
24 Steinfadt et al. 2010. 25 Kilic et al. 2010b . 26 Shporeretal. 2010, 27 Brown et al. 2011. 28 Bassa et al. 2003 . 
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